Enhancement of surface activity in CO oxidation on Pt(llO) 
through spatiotemporal laser actuation 
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We explore the efltect of spatiotemporally varying substrate temperature profiles on the dynamics 
and resulting reaction rate enhancement for the catalytic oxidation of CO on Pt(llO). The catalytic 
surface is "addressed" by a focused laser beam whose motion is computer-controlled. The averaged 
reaction rate is observed to undergo a characteristic maximum as a function of the speed of this 
moving laser spot. Experiments as well as modelling are used to explore and rationalize the existence 
of such an optimal laser speed. 
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I. INTRODUCTION 

A crucial component of chemical engineering process 
design is to locate optimal operating conditions, which 
yield optimized productivity, selectivity and flexibility 
under safety and environmental constraints. Typically 
one searches for optimal steady state conditions, possi- 
bly stabilized through feedback loops; yet it is also well- 
known that non-steady-state operating conditions may 
lead to better average performance Examples of 

such chemical engineering processes include fast pres- 
sure swing adsorption [2l and reverse flow reactors 
Improvement of the average process performance under 
non-steady-state conditions (e.g. periodically varying re- 
action conditions) can sometimes be rationalized through 
certain resonances between the characteristic operation 
time scale (e.g. the period of forcing) and the intrinsic 
time scales of the system Q . 

Catalytic reaction systems often exhibit rich nonlin- 
ear dynamic behavior, including spatiotemporal patterns 
such as propagating reacting fronts, pulses, rotating spi- 
rals and even chaos [1, H, Q- Since the emergence of 
these spatiotemporal patterns is governed by the system 
intrinsic time and length scales, catalytic systems may be 
promising testing fields for optimal non-steady-state op- 
erating policies 8] . The development of advanced surface 
resolving techniques such as PEEM d, EMSI [Xi\, 
and RAM jT2| has facilitated the observation of catalytic 
spatiotemporal patterns at the micron scale. Related 
advances in spatiotemporally addressing catalytic sur- 
faces at the micron scale, in particular through computer- 
controlled laser beams, have enabled us to spatiotempo- 
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rally modify local reaction conditions in real time: pulses 
and fronts, the building blocks of spatiotemporal pat- 
terns, can be initiated, erased or guided through appro- 
priate local laser actuation in space and time . 

In a previous experimental study of modifying the lo- 
cal catalytic surface activity in CO oxidation on Pt(llO) 
through local laser actuation, it was discovered that opti- 
mal laser scanning policies exist, maximizing overall reac- 
tion (i.e. CO2 production) rate when the system is close 
to the excitable regime [1J| . The reaction conditions used 
in those experiments were selected in a way such that 
the system was initially at a CO-poisoned state (with Pt 
surface dominated by adsorbed CO which leads to low 
CO2 production rate), but could become locally excitable 
upon local laser heating. The actuation of the laser beam 
converts the local Pt surface to an excitable state by in- 
creasing the local temperature and hence accelerating 
the desorption of CO. When the laser beam is kept at 
a fixed position, reacting oxygen pulses are periodically 
initiated around the laser heating center, and propagate 
away with significant higher local reaction rate. When 
the laser beam is scanned at a constant speed, it contin- 
uously locally excites the Pt surface . Depending on 
the intensity of the laser, the spatiotemporal "policy" at 
which it scans the surface, its scanning speed, and the 
characteristic times required for the surface to return - 
after excitation- to the quenched state, different maxima 
in the overall CO2 production rate were observed and 
partially rationalized in Ref. p^ . 

In order to systematically explore optimal operating 
conditions we choose a simple laser scanning policy (con- 
stant speed "dragging" along a one-dimensional path, 
p^ ) and attempt to establish a quantitative connection 
between the laser scanning policy and the corresponding 
overall CO2 production rate. This is attempted through 
modelling of the appropriate catalytic reaction-diffusion 
system with spatiotemporal laser actuation. As a first 
step, we studied the simple case in which the laser was 
scanned in a fixed circular pattern at various speeds. 



2 



through both numerical bifurcation analysis and simula- 
tions. By identifying a certain pulse instability associated 
with the dragging of the laser spot, our computational re- 
sults could predict the existence of a local maximum in 
overall CO2 production rate while varying the laser scan- 
ning speed, in agreement with the experiments. The pa- 
per is organized as follows. We begin with a brief descrip- 
tion of our model of CO oxidation on Pt(llO) with local 
laser dragging. Computational results from numerical 
bifurcation analysis and transient simulations are then 
presented and discussed. We also briefly describe the ex- 
perimental setup and present the experimental results; 
we conclude with a discussion of the relation between ex- 
periments and modeling, and possible extensions of the 
work. 



II. MODELLING 

The experimental system of interest is the UHV oxi- 
dation of CO on Pt(llO); we therefore use the Krischer- 
Eiswirth-Ertl (KEE) reaction-diffusion model for this re- 
action p^ . The surface reaction follows a Langmuir- 
Hinshelwood mechanism: 

C0 + * ^ COad., 
2*+02 ^ 20ad., 
COads + Oads ^ 2*+C02T, 



accompanied byalx2-^lxl phase transition of the 
Pt(llO) surface due to CO adsorption. The equations for 
this kinetic model are 



■ii = KsuPco 



U.9 



klu — k2uv + V • (-D„Vu), 
(1) 
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kouv. 



w = k3{f{u) - w), 



(2) 
(3) 



where u, v and w denote the surface coverage of CO, O, 
and the surface fraction of the 1x1 phase respectively. 
The adsorption rate constants for CO and O2, ku 
and ky respectively, are set to fixed values within the 
temperature range considered in this paper. The rate 
constants fci , k2 and k^ for the desorption, reaction and 
surface phase transition are temperature dependent 
through the Arrhenius formula ki — kf eicp{—Ei/kT); 
T is the temperature of the single crystal. We used the 
parameters for Pt(llO) given in Table II of Ref. [31 as 
follows: 



Us = 1, Vs = 0.8, Su = 1, Si,i = 0.6, = 0.4, 

ku = 3.135 X lO'^ s-i mbar-\ 

ky = 5.858 X lO'^ s-^ mbar-i, 

h = k°exp{-E,/kT), i = l,2, and 3, 

fc? = 2 X 10^6 s-\ El = 38 kcal/mol, 

fc§ = 3 X 10^ s"\ E2 = 10 kcal/mol, 

fcO = 10^ s-\ = 7 kcal/mol. 



We adopted the diffusion coefficients reported in 
Ref. : 

Du = Dlexp{-Eu/kT) 

where, in the [110] direction, — 5x lO^'^ cm^/s, Ey — 
10 kcal/mol and in the [001] direction D° = 7 x IQ-^ 
cm^/s, Eu = 8.9 kcal/mol. The function f{u) has been 
fitted to experimental data to give the rate of surface 
phase transition as a function of u, the coverage of CO, 
as follows: 



/(«) 





M^-1.05tt^+0.3M-0.026 
-0.0135 
1 



for u 0.2, 
for 0.2 < M < 0.5, 
for u ^ 0.5. 



In Eqs. (HI)-® the temperature field is assumed to be 
homogeneous across the catalytic surface and prescribed 
by a single parameter T. When a laser beam is scanned 
across the sample surface, the local temperature of the 
catalytic surface around the laser spot is temporarily in- 
creased, leading to spatiotemporal heterogeneity in the 
surface temperature field. In order to model such a spa- 
tiotemporal heterogenous temperature field, we make the 
following two simplifications [H, [l^ : a) the heat gener- 
ated by the CO oxidation reaction on the Pt(llO) surface 
1 mW/cm^) is neglected compared to the power of the 
laser beam (~ 1 W for a crystal surface of 75 fira in diam- 
eter (size of the laser spot)); b) we consider that the lo- 
cal temperature profile inhomogeneity is established (and 
then vanishes) instantly as the laser beam is applied (and 
then removed) . This is because the time scale associated 
with the local temperature change 10^^ second, es- 
timated from W'^/a, where W is the diameter of the 
laser spot (75 fim) and a is the thermal diffusivity of Pt 
(^ 10^^ cm^/s)) is much smaller than the time scales 
of local adsorbate coverage change (~ 10^^ second and 
~ 1 second, estimated from CO desorption with l/ki and 
CO diffusion with /Dy. Based on these two simpli- 
fications, the spatiotemporal change in the temperature 
field induced by a moving laser beam can be modelled 
as a temperature "bump" with fixed shape co-travelling 
with the laser spot. 
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In the experiments, the laser beam is moved in a cir- 
cular path on the sample surface at constant speed. Be- 
cause the typical diameter of the circular path (~ 860 
^m) is significantly larger than the size of the laser spot 
(~ 75 fim in diameter), it is informative to model the 
laser scanning in one dimension with periodic boundary 
conditions. A computational frame co-travelling with the 
laser spot is introduced. Based on the KEE model and 
the discussion above, we have the following equations for 
CO oxidation in a one-dimensional periodic domain with 
laser scanning in a co-moving frame: 



Ut = kuSuPCO 



Us 



kiU — k2UV+{DuUz)z+CUz, 

(4) 

Vt = kvP02[wSvi+{l~w)Sv2] - — j -k2UV+CVz, 



wt = k3{f{u) -w)+ cwz 



(5) 
(6) 



where c is the speed of the co-moving frame (equal to 
the laser moving speed); z is the spatial coordinate in 
the co-moving frame with z = x — ct; and the remain- 
ing symbols are the same as in Eqs. ([H)-©, except fci, 
/c2 and ks , which are now spatially dependent because of 
the nonuniform temperature field T(z) due to local laser 
heating. In the computations the local T(z) heterogene- 
ity due to the laser spot is approximated by a Gaussian 
profile. 
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FIG. 1: (Color online) Enhancement in the CO2 production 
rate as a function of the laser scanning speed for isotropic CO 
diffusion. The steady coverage profile is computed in a frame 
co-travelling with the laser spot. Stable (unstable) steady 
coverage profiles are marked in solid (dashed) lines. Repre- 
sentative coverage profiles are shown in the insets, where the 
blue and green curves represent the coverage of CO and the 
fraction of the 1x1 phase respectively. The coverage of CO 
and the fraction of the 1x1 phase at the center of the laser 
spot are marked by red crosses. Parameters used: domain 
length L — 800, space discretization Aa; — 0.25 with periodic 
boundary conditions. The local temperature field around the 
laser spot is approximated by a Gaussian profile (a — 2), 
with a maximum temperature increase AT = 5K at the 
laser spot center x = 700. Reaction conditions: T — 541K, 
PCD = 4.68 X lO"'^ mbar, poa = 1-33 x 10"* mbar. The unit 
dimensionless time and length correspond to real values of 1 
s and ~ 5.6 ^m respectively. The unit dimensionless speed 
thus corresponds to a real speed of ~ 5.6 nm/s. 



III. COMPUTATIONAL RESULTS 

For simplicity, we first consider the case where the dif- 
fusion of CO in our one-dimensional domain is isotropic 
(i.e. constant £)„). With this assumption, the coverage 
profile of the adsorbates in a static observation frame 
eventually becomes periodic in both space and time, 
when a laser beam is periodically scanned across the do- 
main at a fixed speed c. This oscillatory solution (in 
the static observing frame) corresponds to a steady state 
solution in the co-moving frame (i.e. of Eqs. (H])-®). 
A Newton-Raphson iteration combined with pscudoar- 
clength continuation is used on Eqs. ([I])-® to compute 
and continue the corresponding steady coverage profiles 
as the laser dragging speed, c, is varied (that is also the 
speed of the resulting dragged pulse). The CO2 produc- 
tion rate is then averaged over the entire domain for each 
steady coverage profile and plotted against c, as shown 
in Fig. [1] The reported enhancement is in comparison 
to the reaction rate at the quenched state (without laser 
illumination). Note that the steady profiles of adsorbate 
coverages computed through pseudoarclength continua- 
tion may be unstable (the dashed line in Fig. [T]) , which 
can lead to an estimate of the average CO2 production 



rate very different from direct simulations (through time 
integration) at fixed c. We will discuss this in detail at 
the end of this section. 

In Fig. [U when the laser dragging speed is increased 
from zero to approximately 23 (at point 2), the steady 
profiles of the adsorbate coverage retain a pulse-like ap- 
pearance. Such pulse-like profiles gradually start to "lag 
behind" the laser spot, as the laser dragging speed in- 
creases, which leads to a monotonic increase in the overall 
CO2 production rate to 20% compared to the quenched 
state. Following the solution branch from point 2 to point 
4 in the Figure, the pulse-like structure gradually dis- 
solves and a "plateau" in the CO coverage profile starts 
to develop between the center of the laser spot and the 
lagging pulse-like structure. The "pulse" now appears to 
consist of two distinct fronts - a leading and a trailing 
one. After point 4, the plateau also starts to shrink and 
finally disappears around c = 24. Such a process occurs 
in a very narrow range of c (approximately between 23 to 
24) and causes a very sharp increase in the overall CO2 
production rate enhancement (up to 70%) followed by an 
immediate sharp decrease (down to approximately 5%), 
which is apparently associated with the development and 
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FIG. 2: (Color online) Phase portraits (circles) for steady cov- 
erage profiles at different c values (points f to 6 as in Fig. [1]), 
see text. 



disappearance of a wide plateau in the CO coverage pro- 
file. When the laser scanning speed is increased above 
24, the rate enhancement quickly drops to zero. 

The emergence of the pulse-like structure and the 
plateau in CO coverage profile is thus seen to be inti- 
mately related to the rate enhancement. In order to bet- 
ter understand the nature of this relation, we explore the 
connection between the steady spatial coverage profiles 
and the dynamics of the KEE ODE model, as shown in 
Fig.O It has been argued |2Q,] that, since the oxygen ad- 
sorption is much faster than CO desorption, the temporal 
variation of oxygen coverage can be adiabatically elimi- 
nated for simplification, leading to a reduced two- variable 
model (for u and w). The corresponding nuUclines of this 
two-dimensional model are plotted in Fig. [2] (solid lines). 
Representative dragged pulse coverage profiles in the co- 
moving frame (e.g. points 1 to 6 in Fig. [1]) correspond to 
limit cycles (periodic trajectories) of the ODE system (in 
z) obtained by setting the left side of Eqs. (l4])-([6]) to zero. 
These u and w profiles are superposed on the simplified 
two-equation model nuUclines in Fig. [2] (circles) . 

The trajectories corresponding to points 1 to 3 "reach 
down" to the lower branch of the du/dt = nuUcline, 
confirming the pulse-like nature of the pattern that re- 
sults from laser-dragging. The CO plateau structure, 
clearly seen in the spatial coverage profiles correspond- 
ing to points 3 to 5, are associated with the profile re- 
maining close to the middle branch of the du/dt — 
nuUcline. This "overlap" with the middle branch shrinks 
as the laser scanning speed c is increased. In dynam- 
ical systems terminology, this plateay structure is thus 



associated with a "canard" [2l|. For point 6, the trajec- 
tory does not attain the middle branch of the du/dt = 
nuUcline; it aligns with the dw/dt = nuUcline. 

Based on the discussion above, the relationship be- 
tween the enhancement in CO2 production rate and the 
laser dragging speed c , as shown in Fig. [1] can be par- 
tially rationalized as follows: When c is increased from 
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FIG. 3: (Color online) A blow-up picture of the unstable so- 
lution branch (dashed line) in Fig. [1] Bifurcation points are 
marked by red dots. Following the solution branch from left 
to right, representative eigenvalue spectra (marked by dots) 
corresponding to the labelled points are plotted in the insets 
below. 



zero up to values in the neighborhood of point 2 ('^ 23.3), 
the laser spot locally excites the catalytic surface, creates 
and "drags along" a reacting pulse (e.g. point 1). Such 
a pulse remains attached to the laser spot; it starts be- 
coming more spatially extended as c is increased, leading 
to a monotonic increase in production rate (appreciable 
adsorbed CO and oxygen coexisting over longer spatial 
extents). When c becomes too fast for the laser actua- 
tion to excite the local catalytic surface, the brief tem- 
perature increase effected by the passage of the laser spot 
only leads to a small local variation in the coverage profile 
(e.g. point 6). For values of c between the two limiting 
cases above (e.g. between points 2 and 5 in Fig. [T]), the 
steady coverage profiles computed from pseudoarclength 
continuation become unstable. In this case, direct sim- 
ulations must be used to find what the stable long-term 
behavior is, and thus establish a connection between c 
and the corresponding rate enhancement (see below). 

A blow-up of the unstable solution branch in Fig. [1] is 
shown in Fig. [3] along with the corresponding eigenvalue 
spectra of the system linearization in the neighborhood of 
the bifurcation points. These linearization spectra have 
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been computed in a finite but long domain with periodic 
boundary conditions; computed eigenvalues are marked 
by dots, and the density of these dots is suggestive of 
continuous spectrum components for the infinite domain 
problem. 

As we move from point a to point c, a "parabola" of 
eigenvalues crosses the imaginary axis from left to right. 
For our large but finite periodic problem, the leading 
eigenvalue (the eigenvalue(s) with largest real part on 
this parabola) is a single real one (as in inset (6)). The 
initiation of the crossing of such an "eigenvalue parabola" 
corresponds to a Saddle-Node (SN) bifurcation followed 
by a large number of subsequent Hopf bifurcations. As 
c varies from point c to point d, the leading spectrum 
component appears to split in two. At point e a single 
real eigenvalue (whose path started close to the parabolic 
spectrum component farthest to the left from the imag- 
inary axis) crosses the imaginary axis, corresponding to 
another SN bifurcation. From points / to g, the split 
spectrum components appear merge again into a single 
component, which then crosses the imaginary axis from 
right to left (see point h). The leading eigenvalues of 
this new parabolic spectrum component (the last ones to 
cross to stability) are now a conjugate pair (in contrast to 
point e). The crossing of such a component of the spec- 
trum appears in our finite domain case as (from point g 
to point i) as a chain of Hopf bifurcations. Continuous 
eigenvalue spectrum crossing has been shown to be asso- 
ciated with the emergence of "chemical turbulence" for 
CO oxidation in an excitable medium [22i] . 

For values of c for which the dragged pulse solutions are 
unstable, direct simulation is used to establish the nature 
of the observed dynamics (and to find the correspond- 
ing average reaction rate). Phase diagrams of the long- 
term dynamics for such simulations appear quasiperiodic 
or mildly chaotic, and do not repeat exactly (as shown 
in Fig. U). For sufficiently large domains (as was the 
case here) , the approximate period of the oscillations is 
not affected by the domain length. In this regime, as c 
is increased, the system spends less time on average in 
states with a long, reacting pulse-like tail, which leads 
to a decrease in the overall (space-time averaged) CO2 
production rate. This is consistent with the (spatiotem- 
porally) averaged reaction rate monotonically decreasing 
as c varies from points e through h in Fig. [3] 

Anisotropic CO Diffusion. To incorporate the 
effect of diffusion anisotropy for the adsorbed CO on 
Pt(llO), we approximate the anisotropic diffusion coef- 
ficient with a spatially periodic function, I?„ = {D^f + 
Dus)/2 + {Duf — Dus)/2 X sin{4:TTx / L) , where Duf and 
Dus correspond to the coefficients of CO diffusion on 
Pt(llO) in the fast ([110]) and in the slow ([001]) direc- 
tions respectively. The interaction between laser scan- 
ning and strong anisotropic CO diffusion is illustrated in 
Fig. [5l which shows that the reacting tails only exist in 
regions where the local CO diffusion is slow. To visually 
enhance this phenomenon, and also to show how differ- 
ent the instantaneous pulse profiles can be as the laser 



(a) (b) 




Re 



(c) (d) 




FIG. 4: (Color online) Space-time plots and phase portraits 
of CO coverage dynamics for different laser scanning speeds 
in the unstable regime, (a) and (c), space-time plots of CO 
coverage for c — 23.5 and 23.6 respectively, (b) and (d), the 
corresponding phase portraits of the first fast Fourier trans- 
form component of the CO spatial coverage profile. When c 
is increased the system spends less time on average at a state 
with a long reacting tail; thus the overall CO2 production rate 
is less in (c) as compared with (a). The time interval between 
dashed lines in (a) and (c) corresponds to L/c, that is, one 
"rotation" of the laser beam along the periodic domain. 
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FIG. 5: (Color online) Laser scanning in the presence of 
strong anisotropic CO diffusion, (a) space-time plot of CO 
coverage in a frame co-taveling with the laser spot; (b) repre- 
sentative instantaneous snapshots of the CO coverage profile. 
The anisotropic diffusion of CO is approximated by a peri- 
odic diffusion coefficient with a period of L/2 (dotted line in 
(b)), Du = {Duf ^ D„.)/2 + (D„/ - Dus)/2 x sin(A-Kx/L), 
where D^f and Dus correspond to the diffusion coefficient val- 
ues along the fast and slow directions on Pt(llO) respectively. 
Duf /Dus = 3, c = 23.3. 

rotates on the surface Fig. [5] employs diffusion coefficient 
values six times as large as the ones reported in fl7j (but 
with the correct anisotropy ratio). Resetting D^f and 
Dus to their experimentally reported values [l7| j time in- 
tegrations of Eqs. (HI to ([6]) show periodic or apparently 
quasiperiodic transient behavior (Fig. [7]). In both cases 
we average the reaction rate over both the domain length 
and a long enough time interval to compute the corre- 
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rate is computed based on the spatiotemporal coverage pro- 
file, averaged over tlie entire domain and a sufficiently long 

time interval. D„/ and D^s are estimated from experimental piG. 7: (Color online) Phase portraits for anisotropic CO 
data [17l|; Duf/Dus ~ 3. diffusion at different laser scanning speeds. The local CO 

diffusion coefficient at a; = 700 is plotted on the right of plot 
(a), (a) c = 23.3, (b) c = 23.5, (c) c = 23.6, (d) c = 25. 



spending average enhancement in CO2 production rate 
shown in Fig. [5] When c is varied between and approxi- 
mately 23, the dynamics in the co- moving frame converge 
to a stable oscillatory state with a period of L/ (2c). Dur- 
ing the oscillations, the length of the reacting tail only 
slightly varies as the local CO diffusion coefficient modu- 
lates (not shown). Similar behavior has been observed for 
larger c (i.e. > 24) , where the width of local coverage pro- 
file around the laser spot slightly modulates. Compared 
with the isotropic case, the averaged rate enhancement 
shows the same qualitative behavior for the two limit- 
ing cases: for smaller c the rate enhancement gradually 
increases as c is increased; for large c the rate enhance- 
ment suddenly decreases as c is increased. When c is in 
the intermediate range 23 ~ 24 (associated with insta- 
bility in the isotropic case), the system can show distinct 
quasiperiodic transient behavior: for c close to the lower 
boundary of this interval the system exhibits oscillations 
similar to the corresponding isotropic case (Fig. El^b)); 
for c close to the upper boundary, "skipping" oscillations 
can be observed, as shown in Fig. By comparing 

the position of reacting tails in Fig. [7fc) with the corre- 
sponding spatial profile of Z3„, we find that the reacting 
tails occur only along the slow diffusion directions (when 
the local Du is small). This suggests that the local peri- 
odic change in CO diffusion coefficient as the laser spot 
moves along a ring can have a strong impact on the lo- 
cal dynamics; it can, for example temporarily eliminate 
a reacting tail and then reestablish it (see the small time 
gap in Fig. [7{c) when the reacting tail disappears). As 
in the isotropic case, the averaged rate enhancement in 
this range of c monotonically decreases as c is increased. 



IV. EXPERIMENTS 

The simplest motion of a temperature heterogeneity 
on the platinum surface would be a linear one, on an in- 
finitely extended catalyst. This cannot be realized in the 
experiment. The easiest consistent approach is to use 
a protocol where the laser spot moves on a circle with 
large circumference. This is experimentally feasible and 
makes a comparison with theory possible. Circular move- 
ment was tentatively explored in a previous study p^ : 
however, the resulting surface dynamics and changes in 
reaction rate were not analyzed in detail. Here, we show 
new experimental results, focussing on reaction enhance- 
ment and the influence of the laser spot on the reaction 
patterns; we also compare to the theoretical results pre- 
sented in the first part of the paper. 

Experiments are performed using a Pt(llO) single crys- 
tal with 10 mm diameter located in an ultra high vacuum 
(UHV) chamber. For accurate control of the reaction 
parameters (temperature, CO and oxygen partial pres- 
sures), the chamber is equipped with a computer con- 
trolled gas dosing system, that allows for adjusting the 
partial pressures of carbon monoxide and oxygen. The 
platinum sample is heated from the back with a halo- 
gen lamp. For sample preparation and characterization, 
the chamber is equipped with Ar-ion sputtering and Low- 
Energy Electron Diffraction (LEED). A tube (^38 mm in 
diameter) with a cone-like ending (~4 mm in diameter) is 
placed about 2 mm away from the sample surface. It con- 
nects the reaction chamber with a differentially pumped 
Quadrupole Mass Spectrometer (QMS) used for analysis 
of the reaction products. This way, a large fraction of 
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FIG. 8: Schematic view of the experimental setup, not to 
scale. See the text for details. Reproduced from p^ . 



the carbon dioxide produced during the reaction can be 
detected, providing an estimate of the overall surface re- 
action rate. The carbon monoxide signal detected by the 
QMS was stored on a computer. 

Concentration patterns of CO and oxygen on the sam- 
ple surface are imaged at a rate of 25 images per second 
using EUipsomicroscopy for Surface Imaging (EMSI) [1^ . 
They are recorded with a video CCD camera. The back- 
ground of the video signal is subtracted and the signal 
is contrast-enhanced in real time using a downstream 
Argus 20 image processing unit from Hamamatsu and 
stored on a video recorder. The surface reaction can be 
locally manipulated by a Argon ion laser beam which 
is focussed onto the platinum crystal. The spot size 
of the laser is approximately 75 ^m in diameter. Two 
computer-controlled mirrors allow for controlled motion 
of the laser spot on the sample surface. Special motion 
protocols (like linear or circular trajectories of the laser 
spot) can be programmed using Lab View. The partial 
pressures of oxygen and CO are also controlled through 
this software. The experimental setup is sketched in 
Fig. El 

At the beginning of each experiment, the laser beam 
is focussed to a fixed position on the sample surface for 
several minutes, until the temperature of the crystal equi- 
librates. This is necessary because the laser spot not only 
produces a strongly located temperature increase but ad- 
ditionally slightly increases the mean crystal temperature 
(by a few K). The laser beam is then scanned across the 
sample surface with a constant speed for 40 seconds, fol- 
lowing a fixed circular path with a circumference of ap- 
proximately 2.7 mm. This is illustrated in Fig. [9l The 
scanning of the laser beam is stopped for a minute every 
time before a different scanning speed for the laser beam 
is applied; at each speed several measurements are per- 
formed and their results in terms of reaction rate increase 
are averaged. 
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FIG. 9: Circular motion of the laser spot across the platinum 
catalyst, (a)-(c) raw data, (d)-(f) images after additional 
background subtraction, contrast enhancement and cropping. 
The time difference between the snapshots is 0.96 s. The crys- 
tal orientation is indicated in (a). Diffusion along the [110] 
axis is faster (by an estimated factor of 2-3) than in the [001] 
direction. The laser spot is visible as a bright white area in 
the images. Due to the laser motion, a "tail" develops be- 
hind the spot. The reaction parameters are PLaaer=640 mW, 
P02=3.00 X 10"* mbar, pco=8.43 x 10"* mbar, and T = 
513 K. 



In EMSI images, dark surface area is associated with 
more reactive, predominantly oxygen covered surface. 
When the laser spot moves, a dark tail can be observed: 
the laser spot removes some adsorbed CO in its path. 
Since CO poisons the catalytic surface, its removal in- 
creases the catalytic activity. In addition to the one- 
dimensional theoretical observations in the first part of 
this article, we see here the full two-dimensional struc- 
ture of the laser induced reactive surface area. At a cer- 
tain distance from the laser spot, the tail divides, and a 
"swallowtail" composed of two narrow "tongues" can be 
seen. FigurefTUTa) shows the CO2 signal measured by the 
QMS. During the periods of laser spot motion (shaded 
regions) the CO2 signal is increased. The delay between 
stopping the laser spot motion and associated decrease 
of the CO2 production rate is due to the slow time scale 
of the differential pumping of the QMS. 

We measure the increase of the carbon dioxide sig- 
nal as a function of the laser spot velocity. The results 
are shown in Fig. [TUTb). We observe that the average 
reaction rate increases for velocity values up to about 
2.8 mm/s. For higher velocities, the enhancement falls 
precipitously until the laser spot moves so fast that no 
enhancement can be measured. When all other reaction 
parameters are kept fixed, there exists a certain veloc- 
ity of the laser spot at which the enhancement is at it's 
maximum. We call this velocity value the "optimal ve- 
locity" . The experimental curve agrees qualitatively with 
theory: A slow increase of the reaction rate is followed 
by a fast decrease after the optimal velocity has been ex- 
ceeded. This has been observed in our computations both 
for isotropic and for anisotropic diffusion coefficients; in 
the case of anisotropic diffusion the dragged pulse profile 
varies along the laser path because the diffusion coeffi- 
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FIG. 11: Reaction rate enhancement at different experimen- 
tal conditions: Each enhancement curves is normahzed with 
respect to its optimum velocity and arbitrarily shifted with 
respect to the y axis for better comparison. The reaction 
conditions were PLaser=640 mW, po2=3.00 x 10~* mbar, 
Pco=8.60 X 10~^ mbar, and T'=512.5 K for experiments A 
and B, and pco=8.43 x 10~^ mbar, and r=513 K for exper- 
iments C and D. 



FIG. 10: Enhancement of CO2 production, (a) QMS sig- 
nal; the laser spot is only moving during the shaded time 
intervals with a velocity of 0.91 mm/s. (b) Average en- 
hancement for different velocities of the laser spot. The reac- 
tion parameters are PLaser=640 mW, po2=3.00 x 10^* mbar, 
pco=8.60 x 10"^ mbar, and T = 512.5 K. 



cient also varies. The variation in the length of the pulse 
tail in the three snapshots of Fig. [9] is the result of the 
varying orientation of laser motion with respect to the 
crystal axes (and the associated CO diffusion coefficient 
variation) . 

Figure [11] shows the results of experiments at differ- 
ent reaction parameters, normalized with respect to the 
optimal velocity of each experiment. The qualitative be- 
havior of the reaction enhancement appears similar in 
all cases shown. The experimental curves exhibit a clear 
maximum in the overall rate enhancement as the laser 
scanning speed is varied, in agreement with the compu- 
tational results shown above. It is interesting to observe 
certain secondary reaction extrema "on the way" to the 
overall maximum as the laser spot velocity increases, e.g. 
at V = 1.35 mm/s in Fig. llOf b). These were not observed 
in our one-dimensional modeling, and they most proba- 
bly should be attributed to the two-dimensional nature 
of the system. 

As the laser moves on the catalytic surface, it locally 
shifts the system into an excitable state. This results 
in the development of a reactive "tail" behind the laser 
spot. The length and the width of this tail vary with 
the spot velocity as clearly seen in Fig.[T21 Qualitatively, 
the variation of the average tail length with the laser 




FIG. 12: Length of the dark reactive tail at different velocities, 
(a) w=0.94 mm/s, ^=0.23 mm, (b) mm/s, Z=0.38 mm, 

(c) v=1.8 mm/s, /=0.53 mm, (d) v—2.6 mm/s, 1=0.79 mm, 
parameters like in Fig. [9] 



dragging speed mimicks the variation of the reaction rate 
enhancement both in modeling and in experiments: the 
tail initially becomes longer and then precipitously dis- 
appears as the speed is increased. However, due to the 
development of two "tongues" in the back of the tail, 
a clear definition of the tail length in the experiment 
is somewhat difficult. We (rather arbitrarily) chose the 
distance between the laser spot and the point where the 
tail visibly splits as a measure of the tail length. With 
increasing laser speed the tail becomes longer and nar- 
rower. In Fig. [13 (d) the length of the tail has increased 
further, yet the reaction enhancement is almost negligi- 
ble. This appears counterintuitive, since in the theoreti- 
cal part of this article we saw a clear correlation between 
tail length and reaction rate enhancement. However, the 
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one-dimensional modelling ignores the width (and over- 
all two-dimensional structure) of the tail; this makes di- 
rect comparison between the two-dimensional experiment 
and the one-dimensional theory difficult. We expect to 
report on a more quantitative comparison, based on full 
two-dimensional simulations, in a future publication. For 
very high laser speeds in the experiment -as well as in the 
ID model- the reactive tail completely vanishes, since the 
thermal energy deposited per time and surface area in 
the catalytic surface no longer provides a sufficient local 
temperature increase to excite the system. 

V. SUMMARY AND CONCLUSIONS 

Following our initial experimental study on 
improving catalytic surface activity through time- 
dependent operating policies, we performed here a sys- 
tematic exploration of the effect laser spot motion on 
the average reaction rate of CO oxidation on a Pt(llO) 
surface. Quantitative relations between overall rate 
enhancement and the laser scanning speed, for both 
isotropic and anisotropic CO diffusion, have been ob- 
tained through numerical bifurcation analysis and direct 
simulations. A characteristic local maximum in the rate 
enhancement (gradual increase followed by precipitous 
drop) has been observed in both cases. The local maxi- 
mum for isotropic CO diffusion was associated with the 
development of a pulse instability caused by laser drag- 
ging and apparently involving the dragged pulse continu- 



ous spectrum. When anisotropic CO diffusion is included 
in the model, and for a wide range of laser scanning 
speeds, the system exhibits comparable averaged rate en- 
hancement response; yet the anisotropy has a significant 
qualitative effect on the detailed dynamics (including 
quasiperiodicity and even mild spatiotemporal chaos). 
The computational results arc supported by experimen- 
tal observations showing a qualitatively comparable in- 
crease of averaged reaction rate with laser scanning ve- 
locity. Optimal velocities were obtained experimentally 
for a range of different reaction conditions. Experimen- 
tal observations revealed an interesting two-dimensional 
structure of the dragged pulse "reactive tail" ; this clearly 
suggests that full two-dimensional modeling, taking into 
account diffusion tranverse to the pulse motion, is re- 
quired for a more quantitative comparison of theory with 
experiments. We are currently working towards such a 
fully two-dimensional model study. The laser motion ex- 
plored in this paper was, in some sense, the simplest spa- 
tiotemporal one: effectively one-dimensional with con- 
stant speed. Spatiotemporally more complex motions, 
such as those presented in possibly also containing 
a stochastic component [24|, should provide a rich field 
of study. 
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